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ABSTRACT 

Elastic-plastic finite-element analyses were used to calculate 
stresses and displacements around a crack in finite-thickness plates for 
an elastic-perfectly plastic material. Middle- and edge-crack specimens 
were analyzed under tension and bending loads. Specimens were 1.25 to 20 
mm thick with various widths and crack lengths. A global constraint 
factor Qg, an averaged-normal -stress-to-flow-stress ratio over the plastic 
region, was defined to simulate three-dimensional (3D) effects in two- 
dimensional (2D) models. For crack lengths and uncracked ligament lengths 
greater than four times the thickness, the global constraint factor was 
found to be nearly a unique function of a normalized stress-intensity 
factor (related to plastic-zone-size-to-thickness ratio) from small- to 
large-scale yielding conditions for various specimen types and thickness. 
For crack-length-to-thickness ratios less than four, the global constraint 
factor was specimen type, crack length and thickness dependent. Using a 
2D strip-yield model and the global constraint factors, plastic-zone sizes 
and crack-tip displacements agreed reasonably well with the 3D analyses. 
For a thin sheet aluminum alloy, the critical crack-tip-opening angle 
during stable tearing was found to be independent of specimen type and 
crack length for crack-length-to-thickness ratios greater than 4. 

INTRODUCTION 

The importance of constraint effects in the failure analysis of 
cracked bodies has long been recognized by many investigators. Strain 
gradients that develop around a crack front cause the deformation in the 



local region to be constrained by the surrounding material. This 
constraint produces multiaxial stress states that complicate stress 
analyses and influence fatigue crack growth and fracture behavior. The 
level of constraint depends upon the crack configuration and crack 
location relative to external boundaries, the material thickness, the type 
and magnitude of applied loading, and the material stress-strain 
properties. In the last few years, a concerted effort by the fracture 
community (see refs. 1 and 2 ) has been undertaken to quantify the 
influence of constraint on fracture. To evaluate various constraint 
parameters, two- and three-dimensional (2D and 3D) stress analyses have 
been used to determine stress and deformation states for cracked bodies. 

In the past, Levy, Marcal , and Rice [3], Kiefer and Hilton [4], de 
Lorenzi and Shih [5], and Moyer and Liebowitz [6] used 3D, small -strain, 
finite-element analyses to study cracked finite-thickness plates. Sih and 
Chen [7] used small -strain plasticity theory but departed from classical 
approaches by varying the tensile flow properties through the thickness. 
Brocks et.al [5] and Horn and McMeeking [9] used large-strain, finite- 
element analyses to study the local stresses and deformations around a 
crack in finite-thickness plates. Similar analyses were also carried out 
by Wilkins [ 10 ] and Ayres [ 11 ] using finite-difference methods and by 
Malik and Fu [ 12 ] with the method of lines. These studies obtained 
plastic-zone shapes and stress distributions around the crack front for 
specific cracked specimens. 

More recently, the T-stress (in-plane tangential stress parallel to 
the crack surfaces [13]), and the Q-stress (stress factor used to 
characterize the difference between actual stresses and the HRR stress 
field [ 14 ]) were used as parameters to predict the effects of in-plane 
constraint on fracture. Analyses by Brocks and Kunecke [J5], and analyses 
and tests by Sommer and Aurich [J6], on compact specimens and on surface- 
cracked plates showed how stable crack growth under monotonic loading 
varied with another constraint parameter, the mean-stress-to-equivalent- 
stress ratio. 

Fatigue-crack growth is also influenced by constraint variations 
around the crack front. A constraint parameter based on the normal - 
stress- to-fl ow-stress ratio (* yy A 0 ) has been used to characterize fatigue 
crack-growth and crack-closure behavior under cyclic loading [ 17 ]. In 
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Newman et.al [J8], a "global" constraint factor a g , the averaged- through - 
the-thickness value of normal stress acting over the plastic region, was 
defined for use in 2D crack analyses to simulate 3D effects. Using the 
global constraint factors, the crack-tip-opening displacements from a 
modified Dugdale model [17] compared reasonably well with 3D finite- 
element results from small- to large-scale yielding conditions [J8]. The 
global constraint factor concept was also successfully applied to predict 
fatigue-crack growth in thin-sheet alloys under aircraft spectrum loading. 

The capability to use fatigue-crack-growth and fracture properties 
from laboratory specimens to predict behavior of structural components 
depends upon the ability to calculate the constraint parameters and to 
determine their effects on critical values of appropriate crack-driving 
parameters (K, J or 5). Unfortunately, parametric studies have not been 
conducted to establish the influence of crack configuration, thickness, 
crack length, applied stress level and material stress-strain properties 
on crack- front stresses, deformations, and constraint for standard 
laboratory specimens and structural components. Because of the complex 
stress fields, such studies must be conducted using 3D numerical stress 
analyses. From these analyses, expressions characterizing constraint need 
to be developed. Given a body with a through-the-thickness crack, 
equations need to be developed which express the crack-front constraint as 
a function of crack configuration, loading, and material properties as 

Constraint = f(c, w, B, S, E, a Q , n) (1) 

These expressions may be used to characterize fatigue-crack growth and 
fracture under various constraint conditions or they may be used in 2D 
models, such as the modified Dugdale model, to correlate and predict 
fatigue-crack growth [17] and fracture [J9]. For use in 2D crack models, 
an averaged-through-the-thickness constraint factor may be required. 

The objective of this paper is to present 3D elastic-plastic, 
finite-element results on the influence of crack length, width, thickness, 
and loading (tension and bending) on a global constraint factor for a 
straight-through crack in an elastic-perfectly plastic material. Standard 
laboratory (middle-crack, double-edge-crack and single-edge-crack) 
specimens were analyzed. Primarily monotonic loading with a stationary 
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crack front was considered, but a thin specimen was also analyzed under 
stable tearing conditions. Comparisons are made between plastic-zone 
sizes and crack-tip-opening displacements calculated from a modified 
Dugdale (strip-yield) model and the 3D analyses to assess the capability 
of the model to simulate 3D effects. Constraint equations are developed 
which relate the global constraint factor to a normalized stress-intensity 
factor for crack-length- and uncracked-ligament-length-to-thickness ratios 
(c/B and b/B, respectively) greater than four. Constraint variations for 
c/B ratios less than four are presented. Measurements made on the 
critical crack-tip-opening angle (CTOA) during stable tearing in a thin- 
sheet aluminum alloy for various specimen types with c/B ratios greater 
than four are also presented. 


NOMENCLATURE 

b Uncracked ligament length, mm 

B Plate thickness, mm 

c Crack length or half-length (see Fig. 1), mm 
d Minimum element size in crack-front region, mm 
E Modulus of elasticity, MPa 

h Plate half-height, mm 

K Mode I (two-dimensional) stress-intensity factor, MPa-/m 
J J contour integral, kN/mm 

M Applied moment, Nm 

n Strain-hardening coefficient 

S Applied stress, MPa 

v Crack-opening displacements along crack surface, mm 
v Q Crack-opening displacement at mid-plane, mm 
w Plate width or half-width (see Fig. 1 ), mm 
x,y,z Cartesian coordinates 

ag Global constraint factor (averaged a yy/° 0 In plastic zone) 
p Constant in constraint equation 

7 Constant in constraint equation 

v Poisson's ratio 

p Plastic-zone size, mm 

p Q Plastic-zone size at mid-plane, mm 

a Normal (opening mode) stress, MPa 
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a Q Uniaxial flow stress, MPa 

6 Crack-tip-opening displacement, mm 

1> c Critical crack-tip-opening angle (CTOA), degrees 

SPECIMENS TYPES AND MATERIAL 

Middle-crack tension, M(T), double-edge-crack tension, DE(T), and 
single-edge-crack bend, SE(B), specimens were analyzed for thicknesses (B) 
ranging from 1.25 to 20 mm with crack-1 ength-to-specimen-width ratio (c/w) 
of 0.3, 0.5 and 0.7 for widths ranging from 10 to 160 mm, see Figure 1. 
Plane-strain conditions were also imposed on some finite-element models to 
show limiting conditions. The material was assumed to be elastic- 
perfectly plastic with properties typical of a high-strength aluminum 
alloy. The modulus of elasticity (E) was 70,000 MPa, the uniaxial flow 
stress (<7 q ) was 500 MPa, and Poisson's ratio (•/) was 0.3. The finite- 
element models were subjected to either monotonic or constant-K loading. 
Under monotonic loading, the crack front was held stationary, but under 
constant-K loading the crack front was grown uniformly through the 
thickness by one element size at a time until the desired crack extension 
was achieved. The stress-intensity factor solutions were obtained from 
reference 20. 


FINITE-ELEMENT ANALYSES 

A three-dimensional (3D), elastic-plastic, finite-element program, 
ZIP3D [21], developed at the NASA Langley Research Center, was used in 
this study. The program uses 8-noded hexahedral elements, the von Mises 
yield criterion, isotropic hardening, small -strain deformation theory, and 
Drucker's associated flow rule. The elastic-plastic analysis is based on 
the initial -stress method and incremental theory of plasticity [22]. 

Crack extension under constant-K loading was modeled by a nodal -release 
procedure [23]. The finite-element formulation and solution method are 
given in reference 21. 

The 3D finite-element model, as shown in Figure 2, contained six 
layers through the half-thickness. The layer adjacent to the mid-plane of 
the specimen (along the x-y plane) is denoted as Layer 1. Thicknesses of 
Layers 1 to 6 were 0.15, 0.125, 0.1, 0.0625, 0.0375 and 0.025 of B, 
respectively. Figure 2 shows the model for a specimen with c/w - 0.5. 
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The width, height, and thickness of the model were scaled to the desired 
values. Region 1 shows the mesh refinement around the crack front. The 
smallest element size, d, in the crack-front region (Region 1) was 
0.0008w. For c/w ratios of 0.3 and 0.7, the model shown in Figure 2 was 
modified by expanding or contracting elements in the regions x < 0.4w and 
x > 0.6w to achieve the required c/w ratio. Thus, the region around the 
crack front was shifted but the element sizes and pattern were unchanged. 
All models had the same number of elements (5,706) and nodes (7,203). A 
uniform displacement was applied across the y - h plane for the tension 
loaded specimens and a bending moment was applied on the bend specimen. 

In the incremental analysis, the displacement (or load) increment selected 
was thirty percent of the displacement (or load) required to yield the 
highest stressed element. Analyses of each specimen type were terminated 
when the plastic-yield region extended across the uncracked ligament. 

RESULTS AND DISCUSSION 

The "global" constraint factor a g is defined as the average normal - 
stress-to-flow-stress ratio for the elements in the plastic zone. In the 
following, the variations of with applied loading for various cracked 
specimens is studied. Comparisons are made between plastic -zone sizes and 
crack-opening displacements calculated from the 3D analyses and from a 2D 
model using the constraint factors for stationary cracks. Constraint 
variations and local normal stresses were also calculated for a crack 
stably tearing in a thin-sheet material. To experimentally study specimen 
type effects on fracture, critical crack-tip-opening angles (CTOA) were 
measured during stable tearing in a thin-sheet aluminum alloy. The 
critical CTOA values are discussed in relation to the constraint factors 
calculated from the 3D analyses. 

Global Constraint Factor 

The global constraint factor was developed for use with existing 2D 
models [17,19] of crack-growth behavior to simulate the effects of 3D 
stress states on plastic-zone sizes and crack-surface displacements. 

These models, which are based on a strip-yield model [24], require that a 
constant stress act in the tensile plastic zone. An average of the normal 
stresses in the plastic zone may be an appropriate value to estimate 
first-order constraint effects. The global constraint factor is 


6 



( 2 ) 


“9 ‘ (1/A T> j/Wl A i 

where is the element area on the y - 0 plane of a yielded element, 
(CTyy/c7 0 )i is the normalized normal stress for the element, and Ay is the 
total area for all elements (N) which have yielded. 

Even for an elastic-perfectly plastic material, the normal stresses 
in the plastic zone are not constant but vary with thickness, loading and 
distance from the crack front. To demonstrate the feasibility of using an 
average constraint factor for the distribution of normal stresses, a crack 
in an infinite plate was analyzed with the strip-yield model. This crack, 
shown in Figure 3, has a nonlinear stress distribution, similar to the 3D 
finite-element results (shown later), acting over the plastic zone. The 
stress distribution was assumed to be 

a yy = *o('*/e) m 

where p* is the plastic-zone size, £ is measured from the physical crack 
tip, and m was selected to be 0.3 to give a nonlinear stress distribution. 
In the strip-yield analysis, the plastic zone was divided into a number of 
segments with a constant stress acting over each segment. The peak stress 
at the crack tip was 3; and the stresses dropped rapidly and approached 
unity at the end of the plastic zone. The average constraint factor from 
equation (2) was 1.4. The plastic-zone size, p, for a * 1.4 was about 80 
percent of p*. (The constraint factor = 1.8 will be discussed later.) 

Figure 4 shows comparisons among the crack-surface displacements 
computed from the strip-yield model for the nonlinear stress distribution 
(solid curve) and those from the model with values of 1.4, 1.8 and 3. 
Results for = 3 underestimated both the crack-surface displacements and 
plastic-zone size (p = 0.2 p*). An value of 1.8 was selected to match 
local crack-surface displacements but the plastic-zone size was still 
underestimated by a factor of 2, as shown in Figure 3. An value of 1.4 
overestimated displacements (dotted curve) but gave a reasonable estimate 
for the plastic-zone size (p = 0.8p*). To fit both plastic-zone sizes and 
crack-tip displacements, a two-parameter model would have been required. 
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However, the objective of this paper is to characterize constraint with 
only one parameter. Thus, an averaged a yy/° 0 was selected. 

Stress States 

From the FEM, the normal stresses, a yy , in the crack-front region 
normalized by flow stress, a Q , are shown in Figures 5 and 6 for thick and 
thin specimens, respectively, subjected to tension and bending loads. The 
M(T) half-width and SE(B) width (w) was 40 mm, and the c/w ratio for each 
was 0.5. Element centroidal stresses along the crack plane are plotted 
for Layer 1 (center) and the two outer layers as a function of x/c. The 
crack front is located at x/c * 1. These results are shown for specific 
values of the normalized stress-intensity factor, K/(a 0 /B), where K is the 
2D stress-intensity factor. (This normalized stress-intensity factor is 
proportional to the square-root of the plastic-zone-size-to-thickness 
ratio.) K was calculated from the well-known formulas for the M(T) and 
SE(B) specimens [20], For the thick specimen (Fig. 5), the normal stress 
was high (about 2.7a Q for the M(T) and 3.5a Q for the SE(B) ) at the crack 
front for all layers except for the surface layer (Layer 6). However, in 
the thin-sheet material (Fig. 6), the stresses at the crack front were 
noticeably lower, indicating a lower local crack-front constraint. The 
central layers of the models showed the highest stresses and the slowest 
drop in stresses away from the crack front. The peak values of normal 
stress at the mid-plane and the stress variations through the thickness 
from the present analyses compared well with those from Kiefer and Hilton 
[4] and Moyer and Liebowitz [5]. 

The horizontal lines in Figures 5 and 6 show the global constraint 
factors, <*g. For the thick material, was lower for the M(T) specimen 
than for the SE(B) specimen (c/B - 1) but the constraint factors were 
nearly the same for the thin material (c/B = 16). As expected, the 
values were lower for the thinner material indicating less constraint than 
the thicker material. 

Specimen Types 

The variation of with stress-intensity factor for a wide range of 
plate thicknesses (B - 1.25 to 20 mm) and three specimen types is shown in 
Figures 7 and 8 for small and large widths, respectively. The c/w ratio 
for all specimens was 0.5. (The stress-intensity factor was normalized by 
a Q /w here because the plane-strain solutions are also shown.) The results 
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from the plane-strain analyses appear to indicate upper limit solutions. 
The dash-dot line at * 1 denotes a lower limit condition. For 
thicknesses less than 5 mm for the 40-mm wide specimens and thicknesses 
less than 20 mm for the 160-mm wide specimens (that is, c/B £ 4), the 
values were nearly independent of specimen type and showed the same 
variation with normalized K. For thicker plates, the M(T) specimen gave 
lower constraint factors than the DE(T) and SE(B) specimens. All curves, 
except for the plane-strain bend specimen, show a loss in global 
constraint as the applied K level (or equivalently the plastic-zone size) 
increases. The plane-strain solutions showed a much slower constraint 
loss with increasing K level. For a given K level, the thinner specimens 
gave the lowest constraint factors. Results for all specimen thicknesses 
tended to approach nearly the same constraint level (about 1.15) as large- 
scale yielding conditions were reached. The rapid loss of constraint at 
K/(* o yw) equal to about 0.9 for the plane-strain M(T) specimen was caused 
by the plastic region extending across the uncracked ligament. The flow 
stress used for these calculations was 500 MPa. Other calculations made 
with a Q equal to 300 and 700 MPa for various thicknesses and the same 
width (not shown) demonstrated that the normalized stress-intensity factor 
K/(a 0 7w) was able to collapse the results onto the same curves. 

Because the global constraint factors were found to be nearly 
independent of specimen type for the thinner materials, as shown in 
Figures 7 and 8, it was suspected that the crack-length-to-thickness ratio 
(c/B) or the uncracked-1 igament-length-to-thickness ratio (b/B) were the 
controlling parameters. The constraint factors for all numerical results 
where c/B and b/B ratios are equal to or greater than 4 are shown in 
Figure 9. The constraint factor is plotted against a normalized stress- 
intensity factor, K/(a Q /B), that is proportional to the square-root of the 
plastic-zone-to-thickness ratio. Plotting in this way, the results for 
all specimen types, thicknesses, widths, and crack-length-to-width ratios 
collapsed onto nearly a single curve. The upper dash-dot lines show the 
limiting results from the plane-strain analyses of the three specimen 
types. The plane-strain solutions are slightly dependent upon specimen 
type at low normalized stress-intensity factors. The lower limit appears 
to be about 1.15 and is nearly independent of specimen type. Results from 
the M(T) specimens tended to suddenly drop towards unity when the yield 
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region extended across the uncracked ligament (see results at K/(ct 0 /B) » 

2, 4 and 8). However, results from the SE(B) specimens showed a very 
slight rise during the development of the plastic hinge (see Fig. 10 at 
K/(ct 0 7B) 53 2 for c/B - 4). The curves are from an equation chosen to fit 
the results for c/B and b/B £ 4. The constraint equation is 

- X 1 ' 5 

a « 1.15 + /? e’ 7N n (4) 

where K n - K/(a Q /B), p « 1.25 and 7 - 0.85 for M(T), p - 1.4 and 7 - 0.95 
for DE(T), and p * 1.55 and 7 * 1.05 for SE(B) specimens. 

For c/B and b/B ratios less than 4, the global constraint factors 
were found to be a function of specimen type, loading, thickness and 
width; and they were not uniquely related to normalized stress-intensity 
factor. Figures 10 and 11 show how the constraint factor varies with 
normalized stress-intensity factor for various widths and c/w ratios, 
respectively, for a thick plate (B * 20 mm). The open and solid symbols 
show the results for the M(T) and SE(B) specimens, respectively. The 
effects of specimen width are shown in Figure 10. The SE(B) specimens 
produced higher constraint factors for smaller widths or lower c/B ratios, 
whereas the M(T) specimens showed a drop in constraint. The curves are 
from equation (4) and are for specimens with c/B and b/B ratios equal to 
or greater than 4. For a given width (w = 40 mm), the effects of c/w 
ratios on constraint are shown in Figure 11. The constraint factors for 
the bend specimen increased with c/w ratios, whereas the results for the 
tension specimen showed very little influence of the c/w ratio. 

Plastic-Zone Sizes 

The variation of the plastic-zone size through the thickness of a 
thin sheet and a thick plate M(T) specimen are shown in Figures 12 and 13, 
respectively. These figures show the normalized plastic-zone size (p/ c) 
along the y - 0 plane plotted against z/B (z - 0 is located at the mid- 
plane of the specimen). The symbols show the finite-element method (FEM) 
calculations for various values of K/(a 0 /w) with w - 40 mm. For the thin- 
sheet material (Fig. 12), the plastic-zone size is largest near the free 
surface at low K values. However, at higher K values the maximum yield 
zone occurs in the middle of the specimen. But further loading causes the 
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yield zone to become nearly uniform through the thickness. These results 
are in contrast to those expected, in that, the maximum plastic-zone size 
would have been expected to be at the free surface. The dashed lines 
denote the plastic-zone sizes calculated from the modified Dugdale model 
(MDM) [17] using the global constraint factors obtained from the 3D FEM. 
The calculated plastic-zone sizes are reasonable average values through 
the thickness. 

Calculated plastic-zone sizes for a thick plate (B - 20 mm) are 
shown in Figure 13 for several values of the normalized stress-intensity 
factor, K/(a 0 7w). Again, the plastic-zone size is largest near the free 
surface at low K values. However, for higher K values the maximum yield 
zone moves inwards towards the specimen mid-plane. The dashed lines 
denote the plastic-zone size calculated from the MDM [17] using the global 
constraint factor (as shown). The calculated plastic-zone sizes are, 
again, reasonable average values through the thickness. 

Figure 14 shows the plastic-zone size at the mid-plane of the M(T) 
specimen normalized by crack length (p Q / c) as a function of the normalized 
stress- intensity factor. The open symbols show the results for a wide 
range of thicknesses. A plane-strain analysis is. also shown (solid 
symbols) for comparison. The plane-strain results show a rapid increase 
in plastic-zone size at a normalized K value of about 0.9. The other 
thicknesses showed similar but less rapid transitions. These transitions 
in plastic-zone size with normalized K for the finite-thickness specimens 
are associated with a rapid loss of out-of-plane constraint. Even the 
plane-strain solution exhibits an in-plane constraint loss when the yield 
region spreads across the uncracked ligament. 

The curves in Figure 14 show the calculated plastic-zone sizes from 
the MDM [17] using values of the global constraint factor, a^, of 1 and 3. 
These calculations tend to bound the plastic-zone sizes on the upper limit 
(<»g « 1) and are reasonable approximations of the lower limit with an 
assumed plane-strain condition (a^ = 3). 

Crack-Tip Displacements 

A comparison of normalized crack-tip-opening displacements (v Q /c) at 
the mid-plane (z/B - 0) of an M(T) and SE ( B) specimen as a function of the 
normalized stress-intensity factor, K/(a Q ,/w), is shown in Figure 15. 
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These results are for the largest width (w - 160 mm) specimen analyzed; 
and the c/B and b/B ratios were equal to or greater than 4. Displacements 
were calculated very near the crack front (second node behind the crack 
front at x/c - 0.9968). The normalized stress- intensity factor correlated 
the 3D FEM displacements for both specimen types even as the plastic zone 
spread across the uncracked ligament. As expected, the thin specimen 
showed larger displacements than the thick specimen at the same normalized 
K level. The solid curves show the calculated crack-tip-opening 
displacements (at x/c * 0.9968) using the MDM [17] with a g values 
calculated from equation (4). The effective flow stress in the MDM was 
a g a o* ^he solid curves agree reasonably well with the finite-element 
results from small- to large-scale yielding conditions. The predicted K 
level at plastic collapse for the thick specimen was about 5 percent 
higher than the 3D FEM results. For comparison, the dotted curves show 
the calculated results from the MDM [17] for various constant values of a 
from 1 to 3. 

Figure 16 shows the crack-tip-opening displacement (v Q /c) at the 
mid-plane of M(T) specimens as a function of the normalized stress- 
intensity factor. The c/B and b/B ratios for these specimens were 0, 1 
and 16 for the plane-strain case, the 20-mm and 1.25-mm thick specimens, 
respectively. The displacements, again, were calculated at the second 
node behind the crack front. The results from the FEM are shown for the 
thinnest and thickest plates analyzed (open symbols). The results for 
plane-strain conditions are shown as the solid symbols. For the same 
normalized stress-intensity factor, displacements were larger for the 
thinner specimen. As expected, the thinnest specimen had a lower 
constraint and, hence, had larger displacements near the crack front as 
compared to the thicker specimen. The plane-strain elastic solution is 
shown by the dash curve. Again, the solid curves show calculated crack- 
tip-opening displacements (at x/c - 0.9968) using the MDM [17] with the a g 
values calculated from the 3D finite-element analyses (as shown by the 
solid curves in Fig. 7). The solid curves agree reasonably well with the 
finite-element results from small- to large-scale yielding conditions. 

Some differences were observed for the thick specimen for K/(<r 0 /w) greater 
than 0.7. The upper and lower dotted curves in Figure 15 show the 
calculations from the MDM for a - 1 and 3, the limiting conditions. At 
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low values of K/(ct 0 7w), the FEM and MDM results approach the plane-strain 
solution. However, at high values of K/(a 0 7w), the FEM results show a 
loss of constraint. The vertical asymtote from the FEM results would have 
been predicted from the MDM using an value of 1.1 for the thin specimen 
and 1.2 for the thick specimen. 

The results shown in Figures 15 and 16 suggest that the global 
constraint factor controls the crack-surface displacements very near the 
crack front; and that the MDM can be used to model three-dimensional 
effects (plastic-zone sizes and crack-surface displacements). 

Stably Tearing Cracks 

The previous analyses were made with stationary cracks. Because the 
normal stresses were extremely high in the crack front region even for 
thin-sheet materials, it was of interest to see how the local constraint 
effects, normal stresses, and crack-surface displacements were influenced 
by stable tearing. Analyses were conducted on a thin sheet M(T) specimen 
where the crack was grown under constant-K conditions. In addition, 
results from previous studies on the influence of specimen type on 
critical CTOA values are presented and discussed in relation to the global 
constraint factors. 

Analyses- -Figure 17 compares the normal stresses for a stationary 
crack and a stably tearing crack in a thin sheet at the same normalized K- 
value. The final crack lengths for both cases is identical. The amount 
of stable crack growth was about equal to the sheet thickness. These 
results show that the normal stresses are not greatly affected by the 
plastic deformations left in the wake of the advancing crack. The normal 
stresses are somewhat lower in the interior region but approach the same 
normal stresses away from the crack front. These results suggest that the 
constraint effects (high normal stresses) are not greatly influenced by 
stable tearing and by the plastic-wake left along the crack surface. 

Experiments - -Photographic techniques have recently been developed to 
measure critical CTOA values during crack initiation and stable tearing in 
thin-sheet materials [ 25 , 26 ], In one method, a high-resolution optical 
microscope was used to record images of the deformed crack surfaces during 
stable tearing. In each frame, the critical value of CTOA was measured at 
several locations (0.3 to 1.3 mm) behind the crack tip. For each frame, 
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the CTOA values were then averaged. In the second method, a digital - 
imaging correlation method [27] was used to record digitized images of 
speckle patterns around the crack-tip location. These two methods gave 
essentially the same CTOA values during stable tearing. A detailed 
description of these methods are given in Dawicke and Sutton [25]. 

Stable crack growth tests were conducted on three specimen types: 
middle-crack tension M(T), compact tension C(T), and three-hole-crack 
tension TH(T) specimens. The standard laboratory specimens, M(T) and 
C(T), were selected to illustrate the influence of specimen type and 
loading on critical CTOA values. The TH(T) specimen was selected to 
measure CTOA in a structurally-configured specimen. This specimen has a 
stress-intensity factor solution similar to a cracked stiffened panel. 

All specimens were made of 2024-T3 aluminum alloy sheet material (B - 2.3 
mm). The critical CTOA values measured on these specimens are shown in 
Figure 18. The critical angle v c is plotted against crack extension, ac 
(visual measurements on surface of specimens). For each increment of 
crack extension, several values of v c were measured and the average is 
plotted. The initial fatigue crack surfaces were flat and showed a small 
amount of tunneling in the interior. During the early stage of stable 
tearing, however, the crack front exhibited severe tunneling in the 
interior (ac * B) and the crack surfaces were still relatively flat. The 
vertical line indicates the approximate location were the crack surfaces 
had completed the transition from flat to shear mode fracture (45 degree 
slant through the thickness). The values measured during crack 
initiation were high but dropped sharply during crack extensions equal to 
about the sheet thickness. Afterwards, the critical angle appeared to 
level off between 5 to 7 degrees. The solid horizontal line is the rb 

c 

value (6 degrees) determined from finite-element analyses [28] to fit the 
maximum failure loads on several M(T) test specimens. Newman et.al [28] 
determined that the high CTOA values found during crack initiation were 
due to the crack-front shape (severe tunnelling in the interior) and the 
large plastic deformations on the specimen surface (recorded with the 
photographic methods). Tests conducted in reference 28 were on relatively 
large specimens; and the c/B and b/B ratios were greater than 4. These 
test results demonstrate that the critical CTOA is independent of specimen 
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type, if the c/B or b/B conditions are met, and further confirm the global 
constraint factor results shown in Figure 9. 

CONCLUSIONS 

Three-dimensional, elastic-plastic (small-strain), finite-element 
analyses were used to study the stresses and deformations around a 
straight-through crack in finite-thickness plates for an elastic-perfectly 
plastic material. Three specimen types: middle-crack tension M(T), 
double-edge-crack tension DE(T), and single-edge-crack bend SE(B) were 
analyzed for thicknesses ranging from 1.25 to 20 mm and widths ranging 
from 10 to 160 mm. Plane-strain solutions were also obtained. For most 
analyses, the crack front was held stationary, but for a thin-sheet 
material the crack front was stably torn. A "global" constraint factor, 
an averaged-through-the-thickness value of normal stress to flow stress 
ratio ( CT yy/ CT 0 ) in ^e P^stic zone (denoted a g ), was defined for use with 
a strip-yield model. Plastic-zone sizes and crack-opening displacements 
near the crack front were compared with predictions from a strip-yield 
model for various thicknesses and applied stress levels. The results of 
this study support the following conclusions: 

(1) The global constraint factor, a^, is nearly a unique function of a 
normalized stress-intensity factor, K/(a 0 /B), for various cracked 
specimen types, widths, and thicknesses when the crack-1 ength-to- 
thickness (c/B) and uncracked-1 igament-length-to-thickness (b/B) 
ratios are greater than 4. An equation relating and K/(a Q ,/B) was 
developed for c/B and b/B ratios greater than or equal to 4. 

(2) For c/B (or b/B) less than 4, the global constraint factor drops for 
M(T) specimens and rises for SE(B) specimens as the c/B (or b/B) ratio 
decreases. 

(3) Using a strip-yield model and the global constraint factors (from 
finite-element analyses or the equation), the plastic-zone sizes and 
crack-tip displacements for the M(T) specimen agreed reasonably well 
with 30 finite-element results for various thicknesses. 


15 



(4) From the finite-element analyses, the maximum plastic-zone size at the 
crack front and along the crack plane occurs In the Interior of the 
specimen (away from the free surface) for both thin and thick 
specimens. 


(5) Local normal stresses, a , are slightly reduced for tearing cracks In 
a thin sheet compared to those calculated for a stationary crack. 

(6) For a thin sheet aluminum alloy, the critical crack-tip-opening angle 
(CTOA) Is constant and is independent of specimen type after a small 
amount of crack extension (c/B and b/B ratios greater than 4). 
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FIG. 1--Specimen types analyzed for finite-thickness plates. 







FIG. 3--Nonl inear- and constant-stress distributions 
for strip-yield model. 



FIG. 4--Comparison of near crack-tip displacements from strip-yield 
model for nonlinear- and constant-stress distributions. 
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FIG. 5- -Normal stresses near the crack front in thick plate 
under tension and bending loads. 



FIG. 6- -Normal stresses near the crack front in thin sheet 
under tension and bending loads. 
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FIG. 7--Global constraint factors for various specimen types 
and thicknesses with small width (w - 40 mm). 



FIG. 8- -Global constraint factors for various specimen types 
and thicknesses with large width (w - 160 mm). 
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FIG. 9--Global constraint factors for various specimen types 
with c/B and b/B ratios greater than 4. 



FIG. 10--Global constraint factors for a thick plate with various 
widths and c/B ratios less than or equal to 4. 
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FIG. 13--Plastic-zone size through the thickness for various 

K/ (cr o yw) levels for an M(T) thick plate with c/w - 0.5 
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FIG. 14--Mid-plane plastic-zone size against K/(a Jyi) for various 
M(T) specimen thicknesses and plane strain with c/w - 0.5. 
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FIG. 15--Mid-plane crack-tip displacements for M(T) and SE(B) 
specimens with c/B and b/B ratios greater than 4. 
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FIG. 16--Mid-plane crack-tip displacements for M(T) specimens 
with various thicknesses and plane strain. 
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FIG. 17--Normal stresses near crack front in a thin M(T) 
specimen with a stationary and tearing crack. 
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FIG. 18--Comparison of measured critical CTOA values for various 
specimen types with c/B ratios greater than 4. 
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